On Transmission Node Structure in Interacting Systems 



J. D. BanE] and C. A. Stafford 
Department of Physics, University of Arizona, 1118 East Fourth Street, Tucson, AZ 85721 

We provide a series of generic results regarding the structure of nodes in the retarded Green's 
function G of an interacting system, as exemplified by the extended Hubbard model. In particular: 
(1) due to an incompatibility between interactions of nearly any form and a precise definition of 
series propagation, degenerate geometric nodes are split or lifted by interactions; (2) degenerate 
nodes generically exist at the boundary between regimes of node splitting and node lifting and, in 
the presence of interactions, they require fine-tuning; (3) degenerate nodes are highly sensitive to 
perturbation and their sensitivity increases with their degeneracy. Moreover, for high degeneracies 
there is a tendency toward lifting rather than splitting. 

We also propose a characterization of the node structure in extended Hubbard models at arbitrary 
filling in terms of either the eigenvalues of G, or equivalently, the roots of a polynomial. This shows 
that "Mott nodes" previously predicted to occur in the transmission spectra of molecular radicals 1 
are fundamentally associated with nodes in the eigenvalues of the retarded Green's function that 
occur in open-shelled systems. This is accompanied by a low-energy two-pole approximation wherein 
each of the eigenvalues of G are mapped onto a Fermi-liquid-like renormalization of the Anderson 
model, for which the exact self-energy is provided. 



I. INTRODUCTION 

Absent interactions, the path integral formalism pro- 
vides an intuitive relationship between the geometry of 
a quantum system and the interference effects exhibited 
by it. However, the nature of this correspondence in in- 
teracting systems is not obvious. 

In this article we consider in particular nodes from de- 
structive interference that arise in the retarded Green's 
function and consequently the transmission function as- 
sociated with extended Hubbard models. Such models^^ 
have been used to, for example, describe transportPWU 
through molecular junctions and mesoscopic systems. 
In this context there is a great deal of experimen- 
tal and theoretical evidence that most nodes present 
within noninteracting models persist in the presence of 
interactionsPEHU 

However, in contrast to this, degenerate nodes present 
in interacting systems were previously observed to 
be split by Coulomb interactions within many-body 
calculations pDES Here we propose that the foregoing ob- 
servation is the result of an inherent incompatibility be- 
tween the presence of interactions and a definition of se- 
ries propagation that can be formulated diagrammati- 
cally. Remarkably, we find violations of this definition 
cause large disruptions of degenerate nodes even when 
there are no interactions between the units arranged in 
series. The reason is fundamental: Amplitudes in inter- 
acting systems are not expressible as a sum of geometric 
Feynman paths. All possible field configuration histories 
contribute coherently to a propagator, and degenerate 
nodes are found to be strikingly sensitive to this (Figure 

We also show for the first time that degenerate nodes 
do exist in interacting systems. However, in light of 
the aforementioned breakdown of series propagation, we 
find that they generically require fine-tuning and that 
their sensitivity to perturbation and tuning increases 
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FIG. 1. (a) Two identical systems, A and B, linked by a 
matrix element V between orbitals 2 and 2 ; . Even with inter- 
actions only within each system, a definition of series prop- 
agation is fundamentally broken by the processes depicted 
in (b). This can cause dramatic changes in the low-energy 
node structure of the propagator, e.g. permitting transport 
between 1 and l' when it is forbidden between 1 and 2. 



with their order. Moreover, we provide the first reports 
that degenerate nodes can be lifted by interactions in- 
stead of split, that they generically lie at the boundary 
between regimes of splitting and lifting, that this split- 
ting or lifting is an ill-conditioned function of parameters 
in the Hamiltonian, and that for high degeneracy there 
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is a tendency toward lifting. 

To provide a simple conceptual framework for under- 
standing the node structure in extended Hubbard mod- 
els, we also develop a simple but remarkably accurate 
functional form for the eigenvalues of the low-energy re- 
tarded Green's function that qualitatively characterizes 
the node structure of the models we consider. It is rem- 
iniscent of two-pole approximations that have been ap- 
plied to the Hubbard model and it is found to be 
equivalent to mapping each eigenvalue of the retarded 
Green's function onto the Green's function of a Fermi- 
liquid-like^ renormalization of the Anderson modelP^ 

Using this, we show that so-called "Mott nodes", re- 
ported previously in theoretical predictions of molecular 
transmission spectra^ and in the Hubbard model, 3 ^ are 
associated with nodes in the eigenvalues of the retarded 
Green's function in extended Hubbard models. We note 
that such nodes also appear at high energies, and iden- 
tify them as being generically due to the interference of 
field configuration histories. In this context, "eigenvalue 
nodes" responsible for Mott nodes are interpreted as be- 
ing due to the destructive interference of particle-like and 
hole-like processes mediated by the same single-particle 
orbital. This refines earlier worlP advancing the inter- 
pretation that the Mott node itself is directly due to 
particle-hole interference, which, while true in the con- 
text previously considered, turns out to be a special case. 

Finally, we also cast the nodes of the retarded Green's 
function as the roots of a polynomial, which provides 
a formal connection between the perturbation of nodes 



and the perturbation of the coefficients of a polynomial. 
While this formalism gives the nodes of the Green's func- 
tion exactly at all energies, we find it practical to develop 
a reduced-order polynomial that characterizes low-energy 
node structure. In this context, split and lifted degener- 
ate nodes correspond respectively to the real and complex 
roots of this polynomial. 

The organization of this article is as follows: In Sec- 
tion [IT] we develop a definition of series propagation and 
remark briefly upon the connection between propagation 
in an isolated system and transmission through it. In 



Section III we provide an overview of the relationship 
between the breakdown of this definition of series propa- 
gation and the disruption of degenerate nodes in the pres- 
ence of Coulomb interactions. In Section [TV] we describe 
a connection between electronic structure and node struc- 
ture in extended Hubbard models that characterizes low- 
energy nodes in such systems at arbitrary filling. This is 
then used to provide insight into a few generic properties 
of the node structure in interacting systems. 

Many of the results that are described in Sections [HI 
and |IV| are based upon a detailed series of case studies 
of nodes in simple extended Hubbard models. For the 
sake of clarity, these do not appear in the main text and 
instead have been organized in Appendix [Aj Similarly, 
the technical details regarding the connection between 
nodes in the transmission function and nodes in the re- 
tarded Green's function appear in Appendix |B| Finally, 
a derivation of the two-pole approximation for the eigen- 
values of the Green's function used in Section [IV| appears 
in Appendix |C) 



II. A DEFINITION OF SERIES PROPAGATION 



Consider the system depicted in Figure [T] consisting of two interacting systems A and B linked by a single matrix- 
element V. For the sake of concreteness, we take an extended Hubbard modePwith spin and particle-hole symmetry 
for the interactions, in which case the Hamiltonian of the composite system A + B is of the form: 

u A+B = - E j^kUv + \ E v$PiPi 

i,3,v i,3 
System A 

v 

System B 

-E^L^+H.cO + ^E^ 5 ^ 

v s 

Coupling 

where pi = ^ a ($ ia di a — 1 and rji = ^ a c\ a Ci a — 1. Here S na and S na respectively create electrons with spin a on the 
nth site of systems A and B . 

The amplitude for a particle added to the ath orbital of system A at time to be observed in the /3th orbital 
of system B at time t is proportional to a Green's function Gp a (t,0) = —ihQ(t)(dp(t)d^ a (0)). We define series 
propagation in this context via the requirement that Gp a be expressible as a coherent sum over Feynman paths, 32 
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each corresponding to a process wherein a particle propagates back and forth through each system in series, eventually 
ending up in the /3th orbital of system B. The most general expression for Gp a under these conditions is: 



Gp a (t,0)=f driG^ a (ri,0)AGf m (t,ri) 
J o 



Direct transmission 



+ fdn f dT 2 f dT 3 Gt(Tl,0)AG* m (T2,T 1 )A*Gt(T3,T 2 )AG% m (t,T 3 ) 

JO J n J To 



Transmission mediated by reflections at times n and T2 



+ ... 



(2) 



where A is an amplitude associated with hopping between A and B, and G A and G B are amplitudes associated with 
propagation within A and B. These do not need to be Green's function associated the isolated systems, though they 
can be. We merely assume all the amplitudes are causal in the sense that G^(£, 0) = for t < 0. 
In this case, time translation invar iance implies: 



/CO 
dnGt^AGf^t-n) 
-co 

+/ 

J — ( 



dn dr 2 dr 3 Gt(Ti)AG* m (T 2 - Ti)A*Gl(T3 - r 2 )AG B Pm {t - r 3 ) 



(3) 



which is a Dyson equation 3 -^ with the self-energy E^(t) = AS(r)Si n Sjm. Applying the convolution theorem then 
yields a geometric series in the energy domain: 



Gpa — G A a L\G B m 



' GnaAGm m A*Gn n AGp m ■ 



1-G B A*G A A 1-G A G B \V\ 2 
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(4) 



The last equality follows from the condition that the hopping amplitude A, which has not been specified till now, 
is equal to V. This ensures that the foregoing definition gives the exact result for noninteracting systems, as the 
self-energy E above is then the tunneling self-energ}El associated with V. For clarity, the energy dependence of the 
amplitudes Gij is left as implicit. 

We note here that, although the intuitive case of Gp a (t,0) = —ihQ(t)(dp(t)d^ a (0)) was considered first to formu- 
late our definition of series propagation, the preceding applies equally to the retarded Green's function from the 
Keldysh formalism: 34 Gp a (t,0) = — iHQ(i)({dp(i), ^(0)}). Transport related quantities are e legantly expressed via 
this quantity, as it combines the amplitudes for both particle-like and hole-like processes ! 11 * 12 * 35 * 

Thus, throughout the remainder of this work we consider only the retarded Green's function, though the preceding 
is not specific to it. The precise relationship between the retarded Green's function of a system and the associated 
elastic transmission function is given in Appendix [B] We also note that when one of the systems is a set of metallic 
electrodes and bare Green's functions are taken for G A,B , the definition presented here reduces to the so-called elastic 
cotunneling approximation that has been studied before P^HSD 

We now consider the relationship between this definition of series propagation and node structure in interacting 
systems. 
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III. TOPOLOGY AND NODE STRUCTURE: 
THE BREAKDOWN OF SERIES PROPAGATION 



Eq. Q implies that a node at energy E in G A a is 
sufficient to cause a node at the same energy in Gp a pro- 
vided G B n is bounded in the vicinity of E. Equivalently, 
within the foregoing formulation of series transmission, 
if a virtual particle with a given energy can not prop- 
agate from a to n, it can not propagate from a to f) 



by way of coupling between n and m. In particular, if 
G A a and G B m both vanish according to a power law, i.e. 
G oc (E — Eq) 71 , then Eq. Q implies the amplitude G a p 
exhibits a degenerate node, vanishing as {E — Eo) 2ri . 

With interactions these properties do not hold because 
the foregoing definition of series propagation, and conse- 
quently Eq. ([2|, breaks down. This is true even absent 
interactions between systems A and B. In the case stud- 
ies presented in section [Al we find that even when there 
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FIG. 2. The element G23 of the retarded Green's function for 
the depicted tight-binding model (£12 = £23 = £34 = t = 3 eV, 
ti4 = A). This system exhibits a degenerate node at E — 
when A = eV (black), but for A = ±0.1 eV the node is 
respectively lifted or split. Although in this case there are no 
interactions, we find that degenerate nodes generically sit at 
the boundary between node splitting and lifting. 




t t = 3eV 



if (n, m) = (1,4) or (n,m) = (4,1) 
otherwise 



FIG. 4. The geometry and extended Hubbard parameters 
associated with the interactiong supernode depicted in Figure 

m 




. . . G A 

. * 1 G B 

• * • G A VG B 

FIG. 3. Diagrams depicting processes that are inconsistent 
with a definition of pure series propagation between system A 
(blue) and system B (green). Diagrams (a) - (c) are consistent 
with a weak definition of series propagation, whereas (d) and 
(e) are not. We find the presence of either (d) or (e) destroys 
degenerate nodes otherwise expected on the basis of series 
propagation. 



are no Coulomb interactions between the systems, and 
even when G A and G B are dressed exactly by the cou- 
pling of A to B, violations the foregoing definition of 
series propagation dramatically alter the structure of de- 
generate nodes at, for example, molecular energy scales. 
Ultimately, the reason is fundamental: Amplitudes in 
interacting systems are not expressible as a sum of geo- 
metric Feynman paths. All possible field configuration 
histories contribute coherently to a propagator. Thus, 
the factorization defined by Eq. Q omits processes like 
those represented by the Feynman diagrams depicted in 
Figure |3j 

Diagrams of class (a) and (b) are present only when 
there are long-ranged Coulomb interactions between A 
and B, and can be interpreted as dressing the Green's 
function in Eq. Diagrams of class (c) can be viewed 
in the same manner but exist even absent intersystem 
interactions. The nonlocal exchange in (d) does require 
long-ranged Coulomb interactions, but can not be treated 
by dressing G A . It is formally equivalent to introducing 
new matrix elements between A and B that alter their 
mutual connectivity. The final class (e) is inconsistent 
with series propagation in the purest sense - it exists 
generically, can not be treated within any formalism that 
involves only single-particle Green's functions for A and 
B, and does not arise from Coulomb interactions between 
the systems. 

In the case studies presented in Appendix [A} we find 
that for two identical systems arranged in series, the 
presence of diagrams of the form (d) or (e) is necessary 
and sufficient for the disruption of degenerate nodes that 
would otherwise be expected on the basis of series prop- 
agation. When present, diagrams of the form (d) usually 
correspond to self-energies that are quantitatively larger 
than those associated with the higher order diagrams (e) . 
However, we find either (d) or (e) is capable of impart- 
ing very large changes to interference effects otherwise 
expected on the basis of series propagation. Under some 
circumstances, we find they compete and lead to novel 
node structure. 

A curious example of this is furnished by an interact- 
ing supernode present in a four-site extended Hubbard 
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tuning associated with the degenerate node in Figure [5j 
the (1,4) element of this Hamiltonian vanishes due to 
cancellation between nonlocal exchange and higher order 
processes. When this occurs H has the same topology as 
consistent with the existence of a degenerate node 
despite the presence of nontrivial Coulomb interactions 
between all the sites. 

As a more general alternative to this topological inter- 
pretation, it is also possible to understand the crossover 
between node lifting and node splitting in terms of elec- 
tronic structure. We pursue this in the next section, ulti- 
mately describing a simple model that can qualitatively 
explain the node structure at arbitrary filling in the ex- 
tended Hubbard systems considered in this work. 

IV. CONNECTING ELECTRONIC 
STRUCTURE WITH NODE STRUCTURE 



FIG. 5. The node locations (a) in a four-site extended Hub- 
bard model presented in section [A] plotted as a function of 
an element U 1 of the interaction matrix. A degenerate node 
exists at E = for a critical value (U 1 « 0.969) but is either 
split or lifted at other values. In the lifting regime, the eigen- 
values of G corresponding to the HOMO and LUMO make the 
dominant contribution everywhere within the HOMO-LUMO 
gap (a-c). In the splitting regime, the dominant contribution 
is from the HOMO - 1 and LUMO + 1 eigenvalues in the 
region bracketed by the nodes (b-d). 



model depicted in Figure [4| which is studied in more de- 
tail in Appendix [Aj In the present context the pertinent 
observation is that, while this system exhibits a higher- 
order node at the particle-hole symmetric point due to 
series propagation when there are no interactions, in the 
presence of Coulomb repulsion this feature is contingent 
upon the fine-tuning of an interaction matrix element. 
When this element deviates from a critical value, the de- 
generate node is rapidly split or lifted (Figure [5] panel 
a). In Appendix [Aj it is shown that the node lifting in 
this case is attributable to diagrams of the form (e) in 
Figure [3] More generally, we find that degenerate nodes 
generically lie at the boundary between the regimes of 
splitting and lifting. 

Similar behavior can be observed in a noninteractnig 
system wherein a degenerate node is either split or lifted 
depending upon the sign of a perturbation to a hopping 
matrix-element (Figure |2j. Indeed, the splitting or lift- 
ing of the degenerate node discussed above can be in- 
terpreted in topological terms by considering an effective 
Hamiltonian defined by H = + Sq 7 . Here is 
the noninteracting part of the extended Hubbard Hamil- 
tonian and Sq 7 i s the Coulomb self-energy evaluated at 
the node location E = 0. In this simple system, such a 
Hamiltonian accurately reproduces the Green's function 
at half-filling over a wide energy range that includes the 
HOMO-LUMO gap. 

Moreover, in Appendix[A]it is shown that at the critical 



The node structure of an interacting system can also 
be understood in terms of electronic structure. To this 
end, it is useful to consider the eigenvalues of the retarded 
Green's function G, which are analogous to transmission 
eigenvalues. In particular, the retarded Green's function 
of a system with no coupling to external degrees of free- 
dom can be expressed exactly as: 

G(£) = y»)(H (5) 

V 

where \v) is an energy-dependent normalized eigenvector 
of G and \ v is the corresponding eigenvalue. In Appendix 
|Cl it is shown that at zero temperature \ v is given exactly 
by an expression of the form: 

K = £ E-AE v + i0+ (6) 

Here the index r] runs over all possible particle- and hole- 
like transitions out of occupied many-body ground states, 
and AErj and Z% are corresponding transition energies 
and spectral weights. 

In principle the number of terms in the sum above 
is huge-at least equal to the number of N ± 1 particle 
many-body states, where N is the filling of the system. 
However, absent orbital degeneracy, at energies near the 
Fermi level each eigenvalue is given approximately by an 
expression of the form: 

K(E) « - _ £ + , Q+ + - - ^ + , Q+ (7) 

where £ p ' h is an effective parameter that gives the energy 
cost of creating a particle or hole with state \v) and Z^ ,h 
is a corresponding spectral weight. This expression is de- 
veloped in Appendix [C] and the justification for the use of 
effective parameters is considered carefully there. Here 
it suffices to note that we find this approximation to be 
remarkably accurate for the systems considered herein. 
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FIG. 6. Panels (a-d) depict the magnitude of the eigenvalues 
of the retarded Green's function for a symmetric two-site ex- 
tended Hubbard model (system A). In (a), (b), and (c) the 
system is filled by two, one, and zero electrons respectively, 
whereas in panel (d) the Green's function is evaluated at the 
Fermi level. Panel (e) shows the the matrix-elements G n m 
of the retarded Green's function, also evaluated at the Fermi 
level. Nodes in G n m correspond to intersections of the eigen- 
values in panel (d) above. 



With long-ranged interactions, e p,h can be approximated 
in terms of noninteracting energies and the charging en- 
ergy of the system; similarly, Z% + « 1. 

Together Eq. ([5| and Eq. ([7]) provide an essentially 
complete qualitative characterization of the low-energy 
node structure of a wide range of extended Hubbard mod- 



els. This is demonstrated by way of example in Figure [6j 
which shows the exact eigenvalues of a symmetric two- 
site extended Hubbard model from Appendix [A) 

The approximation expressed by Eq. ^ corresponds 
primarily to neglecting the correlation-induced reso- 
nances at high energy (those in the shaded region), which 
are inaccessible to a particle propagating at the Fermi 
level. Due to the (anti) symmetry of the eigenvectors cor- 
responding to the bonding and antibonding resonances, 
the nodes in an element of G occur where two of the 
eigenvalues are equal in magnitude. In more general 
cases the correspondence is modestly less intuitive, but 
the framework provided by Eq. ^ still characterizes the 
node structure accessible at the Fermi level, including the 
nodes in G nm that occur at odd fillings on either side of 
the nodes in the eigenvalues. 

In particular, in Appendix[A]it is shown that the nodes 
in G nm in the N = ±1 filling regions occur due to a node 
in one of the eigenvalues \ V1 and that their locations are 
dependent upon the node structure in the half-filling re- 
gion. More generally, such nodes \ v occur in open-shelled 
systems. Like the correlation induced nodes at high en- 
ergies, they are a signature of non- geometric destructive 
interference between field configuration histories. 

More formally, Eq. (|6| implies the energies of nodes in 
G nm are given exactly by: 



E 



A2 



E - AE V + i0+ 



(8) 



where A v = ZJj (n\u) (u\m) and r] ranges over all transi- 
tions between many-body ground states and many-body 
states with an extra particle or hole. Putting the terms 
in this sum above over a common denominator yields the 
expression: 



A£ 7 ) = 



(9) 



where we have neglected infinitesimal imaginary terms. 

The order of this polynomial is huge because it de- 
scribes the location of all nodes, i.e. even those due to 
correlations at very high energies. Applying Eq. ^ ren- 
ders it more intuitively useful, in which case we have: 



E^ h ii^ 







(10) 



where both \i and v range over single-particle eigen- 
vectors. The degree of the polynomial is then at most 
2{N — 1), where TV is the number of single-particle eigen- 
states. If we take into account the considerations in Ap- 
pendix [Cj it can be seen that a polynomial of order TV — 1 
suffices when the system is in a closed-shelled state. In 
practice, only a few of these roots are present at low- 
energies. Put another way, the number of nodes at low 
energy is at most half the number of eigenvalues compet- 
ing in this region. 
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In light of the present discussion, the node structure 
depicted in Figure |5j explained topologically in the pre- 
ceding section, can also be understood either in terms of 
the interference of the eigenvalues of G (panels a-c) or in 
terms of the roots of the polynomial in Eq. (10). In the 



former case the crossover between lifting and splitting 
occurs as the pair of eigenvalues making the dominant 
contribution to G near E = switches from the HOMO 
and LUMO to the HOMO - 1 and the LUMO + 1. In 



the latter view, lifting occurs when Eq. (10) has com- 
plex roots whereas the splitting regime corresponds to 
real roots. 

These equivalent perspectives generalize readily to the 
case of more complex systems, and imply that the node 
structure in similar models satisfies a few generic condi- 
tions: (1) there is a strong tendency toward the preserva- 
tion of the parity (i.e. evenness or oddness) of the num- 
ber of nodes] (2) all degenerate nodes lie at the boundary 
between the regimes of node splitting and node lifting] (3) 
nodes in the eigenvalues of the Green's function occur 
in open-shelled states and at high energies due to non- 
geometric interference between field configuration histo- 
ries and are responsible for the Mott nodes in G nm ] (4) 
the sensitivity of a node to perturbation or the tuning of a 
parameter increases with its degeneracy] (5) the tendency 
of degenerate nodes to lift rather than split increases with 
their degeneracy. 

The last two points can be understood heuristically 
from the geometry of the eigenvalues \ v . Alternatively, 
they can inferred formally based on Eq. (10) and well- 



known worlP2 regarding the response of polynomial roots 
to the perturbation of their coefficients. In particular, the 
degenerate roots of a polynomial are an ill-conditioned 
function of its coefficients. The same work implies that 
even the location of nondegenerate nodes may in prin- 
ciple be an ill-conditioned function of parameters in the 
Hamiltonian, although this is expected to be a rare case. 

These properties are exemplified by the extreme sensi- 
tivity of degenerate geometric nodes to any perturbation 
inconsistent with series propagation. For example, for 
two systems arranged in series that do not interact with 
each other but are themselves interacting, we find in Ap- 
pendix [A] that processes fundamentally inconsistent with 
the notion of series propagation are sufficient to dramati- 
cally destroy degenerate nodes. Remarkably, in the same 
models these small deviations from pure series propaga- 
tion have little effect on other low-energy features. 



V. CONCLUSIONS 

The retarded Green's function of an interacting sys- 
tem has the potential to exhibit rich and nontrivial node 
structure. However, in this work we have made a few 
generic observations: In extended Hubbard models (1) 
there is a strong tendency toward the preservation of 
the parity (i.e. evenness or oddness) of the number of 
nodes; (2) degenerate nodes require fine-tuning in the 



presence of interactions and sit at the boundary between 
regimes of node splitting and node lifting; (3) nodes in 
the eigenvalues of the Green's function occur in open- 
shelled states and at high energies due to non-geometric 
interference between field configuration histories; (4) the 
sensitivity of a degenerate node to the tuning of a param- 
eter increases with its degeneracy; and (5) the tendency 
of degenerate nodes to lift upon perturbation rather than 
split increases with their degeneracy. 

These properties can be understood in terms of elec- 
tronic structure by way of a simple approximation that 
qualitatively explains the node structure of the extended 
Hubbard models we considered, regardless of filling. 
More formally, nodes of the retarded Green's function are 
the roots of a polynomial. In some cases, node structure 
can also be understood in topological terms. However, 
in the presence of interactions, we find that geometric 
degenerate nodes predicted to exist on the basis of series 
propagation are present if and only if a definition of se- 
ries propagation formulated herein is imposed artificially. 
This is true even when there are no interactions between 
the units arranged in series. 

The reason for this is fundamental: Amplitudes in in- 
teracting systems are not expressible as a sum of geomet- 
ric Feynman paths. All possible field configuration histo- 
ries contribute coherently to an interacting propagator. 
Among the low energy interference features in interac- 
tions systems, degenerate nodes appear to be uniquely 
sensitive to this. 
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Appendix A: Case studies: Detailed analyses of the 
node structure in select extended Hubbard models 



In the following we present detailed analyses of the 
node structure in select two and four site extended Hub- 
bard models. We also touch briefly upon the node struc- 
ture in larger systems. The material presented herein 
forms the basis for many of the observations made earlier 
in this work, but is necessarily more detailed and more 
technical than what has come before this point. The 
reader is encouraged to study the portions of interest to 
him or her. 

Our starting point is Figure [7| which depicts a pair of 
two-site systems, A and B, joined by a single hopping 
matrix element V. The interactions are of the extended 
Hubbard form so that the full Hamiltonian is: 

/ H AJrB = Humana dma + - U nrn p n p m (Al) 



where p n = S n d n — 1, which ensures particle- hole 
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FIG. 7. Top (a): Two sites in an extended Hubbard model 
linked by a matrix element t. Bottom (b): Two systems of 
the form in (a) linked by a matrix element V and possibly 
with some Coulomb interactions between them. This is a 
realization of Figure [I] if we identify sites one and two as 
system A and sites three and four as system B. 
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This system has been studied as a model of a cross- 
conjugated molecule elsewhere, 1 but here we do not con- 
cern ourselves with detailed models of electronic struc- 
ture that might obscure conceptual issues. Instead, we 
take a deliberately simple parameterization and tinker 
with it. In particular, we begin with: 



eV, t = 3 eV, 
Up = 4 eV, 



U = 10 eV, U a = 6 eV, 
V = t, U 1 = U a (A2) 



This parameterization is roughly consistent with in- 
teractions that might be expected on the molecular 
scale JSl^MH in particular for systems with a geometry re- 
specting the symmetries of (b) in Figure [7] 



1. Two sites (System A) 

We first consider system A in isolation. The elements 
G nm of the retarded Green's function in this case are 
plotted in Figure [6] panel (e) in the zero-temperature^ 



limit but with the system maintained at a fixed chemi- 
cal potential /i set equal to the energy E. Here, /i sim- 
ply dictates the filling of the system. Panels (a-d) show 
the eigenvalues of the retarded Green's function either 
at fixed filling (a-c) or with the chemical potential equal 
to the energy (d). The shaded blue areas are the regions 
that are inaccessible to particles propagating at the Fermi 
level, i.e. the regions that do not contribute to linear re- 
sponse transport! 35 * 43 * 

We now comment upon several noteworthy features in 
this figure: 

a. The node in G22 at E = eV: The noninter- 
acting Green's function (not shown) exhibits a node at 
half-filling at E = 0. Formally, this is an antiresonance 
that occurs at the on-site potential of site one regardless 
of the matrix element t. Alternatively, it manifests as de- 
structive interference between propagation mediated by 
the bonding and antibonding orbitals, i.e. the eigenstates 
of the Green's function with positive and negative parity 
respectively. 

In panel (e) it is evident that this node persists in the 
interacting Green's function G22- In fact, there is a sim- 
ple symmetry argument that requires this: While in the 
interacting case there are no single-particle orbitals, the 
eigenvectors of G are still dictated entirely by symmetry 
and therefore identical to their noninteracting counter- 
parts. In the present case, this observation is sufficient 
to imply that the node in #22 must also exist in £^22- Al- 
ternatively, this co uld be viewed as a consequence of the 
Luttinger theorem! 1 ! 44 * 45 ! 

b. Singularities in the bonding and antibonding eigen- 
values: In Appendix [Cj the eigenvalues \ v are given 
by the coherent sum of amplitudes associated with pro- 
cesses wherein a particle or hole is added to a many-body 
ground state. Absent interactions, each of these processes 
is simply the addition of a particle or hole to a nonin- 
teracting orbital, and the energy cost for corresponding 
particle-like and hole-like processes differs only in sign. 
Under these circumstances, the resonances for both kinds 
of processes occur at the same energy and there is one 
singularity per eigenvalue. With interactions this breaks 
down for two reasons: (1) When the system is charged 
(here, away from half filling) the energy cost of adding 
a particle to the system is no longer equal in magnitude 
to that of a hole; (2) transitions between ground states 
and correlated excited states lead to narrow resonances 
at high energies. 

Reason (1) gives rise to Coulomb blockade^M^ in the 
context of transport, and occurs generically in presence 
of repulsive long-ranged interactions. It can be examined 
separately from (2) by taking U nm = U, in which case the 
repulsive Coulomb energy depends only on the net charge 
of the system and the eigenvectors are exactly equal 
to the noninteracting ones. In either cases (1) or (2), each 
singularity corresponds to some physical process wherein 
a particle or hole is added to the system in a manner that 
respects the symmetry of the corresponding eigenvector. 
This leads to novel node structure, as described below. 
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c. The node in the bonding eigenvalue at N = 1 
near E = —10 eV: When N = 1 the bonding orbital 
is half-filled and thus accommodates both particle-like 
and hole-like addition. As noted above, in the presence 
of Coulomb interactions the resonances associated with 
these processes occur at distinct energies. Midway be- 
tween them [E ~ — 10 eV) the energy needed to create a 
virtual particle in this orbital is equal to that of a virtual 
hole, and so these processes contribute opposite ampli- 
tudes, completely suppressing propagation mediated by 
this orbital. 

Note that the crucial fact here is not so much that 
the interference is between particle-like and hole-like 
propagation-this is formally responsible for nodes in G 
including the antiresonance in G22~t>ut that the destruc- 
tive interference is between two processes with different 
amplitudes that are mediated by the same orbitalPThis 
is totally alien to noninteracting systems and hence is be- 
yond the scope of geometric Feynman paths, instead aris- 
ing from the interference of field configuration histories. 
Nodes of this form are actually a general phenomenon 
that can be seen to arise between other singularities in the 
eigenvalues of G. Such nodes correspond to the coherent 
interference of histories (processes) involving states with 
the same symmetry. In a noninteracting system, they do 
not exist. 

d. The nodes in G nm at N = 1 near E = —10 eV: 
The aforementioned nodes in the eigenvalues of G are 
associated with nodes that occur in the elements of G 
around the same energy. Related nodes have been re- 
ported previously in the theoretical transmission spec- 
trum of molecular radicals^ and in the Hubbard modelP^ 
In the present context, a node in the bonding eigenvalue 
can be seen around —10 eV in the N — 1 region. At this 
energy the positive parity eigenvector makes no contri- 
bution to the Green's function and G is perfectly anti- 
symmetric, i.e. Gn = — Gi2- This feature is bracketed 
by nodes in G\\ and G\i that occur due to destructive 
interference between propagation through the bonding 
and antibonding resonances, i.e. at the locations where 
the magnitude of the eigenvalues intersect in panel (d). 
Between these nodes, propagation is mediated almost en- 
tirely by the antibonding resonance, despite the proxim- 
ity of the bonding resonances. 

Considering energies from left to right in the figure, 
as the bonding resonance is suppressed all the nodes in 
G nm in the half-filling region reappear to the left of the 
node in the eigenvalue around —10 eV. As the bonding 
resonance opens back up, one node appears for each ele- 
ment of G that does not exhibit a node in the half-filling 
region. Thus the location of a node in G nm at N = 1 
is related to the existence or nonexistence of a node in 
the same element of G nm at half-filling. An exceptional 
case is when the suppressed eigenvector is the only one 
that contributes to a particular element of the Green's 
function; in this situation the node in that element of G 
coincides with the node in the eigenvalue. Till now this 
is the only case that has been studied, but here we point 



out that it is the exception rather than the rule. Overall, 
this behavior generalizes to larger systems, although the 
structure of the eigenvalues there is richer and can be 
complicated by, e.g., avoided crossings. 

We now consider in a similar manner the node struc- 
ture associated with the composite system A + B. 



2. Four sites (System A + B) 

We now consider the composite system depicted in Fig- 
ure [7] (b) wherein system A is coupled to an identical 
system B by a matrix element V, forming a composite 
system A + B. The node spectrum in this case is de- 
picted in Figure [8] in a manner analogous to Figure [6j 
which depicts the node structure of system A alone and 
is discussed at length in the foregoing section. 

Most of the observations made there carry over to the 
case of the four-site system under consideration here. For 
example, the node structure at odd filling is related in the 
same way to nodes in the eigenvalues of the Green's func- 
tion. However, for the remainder of this section we focus 
our attention on the case of half-filling. In particular, the 
node structure of G23 is noteworthy here: 

Without interactions (not shown) the Green's function 
of system A + B is given by Eq. Q. In this case, the 
node at E = in G22 in system A and its counterpart 
in system B give rise to a degenerate node at E = in 
system A + B where G23 oc E 2 . However, as reported 
previously, 26 in the presence of interactions this is not 
the case. Instead, the degenerate node is split into two 
ordinary nodes which appear around ±4 eV in panel (e) 
for the parameterization given by Eq. ( |A2[ ) . Since this is 
contrary to the predictions of Eq. Q, this node splitting 
is inconsistent with the definition of series propagation 
discussed earlier. It is thus attributable to some amalgam 
of the diagrams in Figure [3] 

To shed light on what is happening, we use exact di- 
agonalization to isolate several interesting combinations 
of these diagrams as well as consider a variety of qual- 
itatively inter estin g variations on the parameterization 
given by Eq. (A2). Along the way, we demonstrate re- 



sults highlighted earlier in this work via a careful con- 
sideration of these cases. In particular, we show that the 
node splitting phenomenon just discussed occurs because 
degenerate nodes are extremely sensitive to the break- 
down of series propagation, even when there are no in- 
teractions between systems A and B. 

The cases considered are organized into the panels in 
Figure [9j which we now remark upon individually: 

a. Node structure with no diagrams from Figure [3J* 
By calculating the exact Green's functions of systems 
A and B in isolation and inserting these into Eq. Q, 
we may calculate the Green's function of system A + B 
excluding precisely the diagrams depicted in Figure [3] 
This is equivalent to enforcing the stritest definition of 
series propagation proposed in section[Hj i.e. one wherein 
G A and G B in Eq. Q are not dressed by the coupling 
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^| Inaccessible at E = u 
~^ Accessible at E = u 
[aj [b^ [c7] [d^ Magnitude of eigenvalues of G 
[eT| Magnitude of elements of G 



FIG. 8. Panels (a-d) depict the magnitude of the eigenvalues 
of the Green's function for system A + B in a manner anal- 
ogous to Figure 6] The associated eigenvectors (not shown) 
are now energy-dependent, but can still be classified based on 
symmetry. Panels (e) and (f) show the elements of the inter- 
acting and noninteracting Green's functions G nm - The G23 
element exhibits a pair of nodes that arise from the splitting 
of a degenerate node present in the absence of interactions. 
In the odd-filling regions a node in an eigenvalue of G causes 
clusters of nodes in G n m- 



of system A to system B. 

The Green's function of A + B calculated in this man- 
ner is shown in panel (a), and a degenerate node can be 
seen in G23. This occurs because the elements an d 
G B 3 of the isolated Green's functions of A and B (not 
shown) exhibit nodes at E = 0, and, with series propaga- 
tion enforced, these necessarily give rise to a degenerate 
node in the composite system. 

We now probe the effects of including one or more di- 
agrams that are inconsistent with at least one of the def- 
initions of series propagation considered earlier. 

b. Node structure with diagrams of the form (a) + 
(b): A priori, (a) and (b) in Figure[3]are plausible candi- 
dates for the splitting of the degenerate node in G23, and 
we can include these diagrams selectively via exact diag- 
onalization. For combinations of (a) and (b) this entails 
calculating G A and G B including intersystem Coulomb 
interactions but with V = 0. Inserting these into Eq. Q 
then gives the desired Green's function, which is consis- 
tent with a weak definition of series propagation wherein 
qA,b are dressed to account for correlations caused by 
long-ranged Coulomb interactions between A and B. 

The Green's functions G A,B calculated during this pro- 
cedure (not shown) still exhibit nodes at E = 0. Con- 
sequently, the the Green's function of A + B calculated 
in this manner, which is depicted in panel (a), has a sec- 
ond order degenerate node. Thus, correlations induced 
by long-ranged Coulomb interactions between A and B 
are not responsible for splitting the degenerate node in 
this case. 

c. Node structure with diagrams of the form (c): Di- 
agrams of the form (c) can also be included selectively 
using exact diagonalization. To this end, the full Green's 
function of A + B is first calculated without long-ranged 
Coulomb interactions but with V 7^ 0. The correspond- 
ing Coulomb self-energy is then determined using the 
Dyson equation, and the parts of this associated with 
self-energies for systems A and B are extracted. These 
are used to determine the Green's functions G A,B dressed 
by the diagram shown, which are in turn used in Eq. Q. 

The elements G A 2 an d G B 3 of the dressed Green's func- 
tions of systems A and B that are obtained during this 
procedure (not shown) still exhibit nodes at E = 0. Con- 
sequently, the degenerate node in the Green's function of 
system A + B persists. 

d. Node structure with diagrams of the form (a) - (c) 
Combinations of diagrams (a)-(c) may be selected for via 
exact diagonalization using almost the exact same pro- 
cedure as for diagrams (c) alone. The only change neces- 
sary is to include long-range Coulomb interactions in the 
initial exact calculation of the Green's function of system 
A + B. 

Again, the dressed Green's functions of A and B (not 
shown) possess nodes at E = that cause degenerate 
nodes in the Green's function for system A + B when it 
is calculated using Eq. (@. 

e. Node structure with diagrams of the form (a) + 
(d) (self- consistent Hartree-Fock): To probe the effects 
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[a^| Series 

[c~| Series with tunneling correction 
[eT| Self-consistent Hartree-Fock 
[g^| Hubbard 

Constant self-energy 
[k^| Hubbard with weak intersystem hopping 
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[b^| Series with intersystem interactions 
[d?| Series with exact self-energy for A and B 
[f^| No interactions between A and B 
[h^| Constant interactions 

System B strictly noninteracting 
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FIG. 9. The absolute elements of the retarded Green's function of system A + B at half- filling, calculated via exact diagonal- 
ization. The calculation is performed with various parameterizations or in such a manner as to selectively include the processes 
depicted schematically in Figure [3] In this case, the presence of diagrams of the form (d) or (e) is necessary and sufficient to 
destroy the degenerate node in G23. Other features in the low-energy region accessible to particles at the Fermi level (white) 
show no such sensitivity here. 
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of diagrams of the form (a) and (d), self-consistent 
Hartree-Fock calculations were performed to determine 
the Green's function of system A + B. The degenerate 
node is split, and calculations including only the Hartree 
diagram (not shown) indicate that the Fock diagram is 
responsible. 

Thus nonlocal exchange is sufficient to split the de- 
generate node; however, it is not necessary or even the 
dominant contribution to the node splitting, as the re- 
maining cases demonstrate. 

/. No interactions between systems A and B: Here 
the Green's function of A + B is calculated exactly with 
no interactions between A and B (i.e. XJab — 0)- The 
degenerate node is split by roughly 3 eV and, on the 
basis of the foregoing cases and the fact that there are 
no intersystm interactions whatsoever, we conclude that 
diagrams of the form (e) are responsible. 

g. Hubbard interactions: Here the Green's function 
of A + B is calculated exactly with Hubbard interactions 
(i.e. U nm = S nrn U). Again, a very large splitting is 
present despite no long-ranged interactions whatsoever. 
In light of the other cases, we conclude diagrams of the 
form (e) are responsible. 

h. Constant interactions: As it is mentioned repeat- 
edly in this work, we consider here the case of constant 
interactions, i.e. U nm = U where U is the average of 
the full interaction matrix from Eq. ( A2 ) . Again a large 
splitting is present in G 2 3- 

i. Fixed Coulomb self- energy: As it is relevant else- 
where in this work, we present here the case wherein the 
Coulomb self-energy E c is fixed to its value at E = 0. 
The splitting shown is attributable to the elements of 
the Coulomb self-energy that connect systems A and B. 
These are precisely the elements that arise from diagrams 
of the form (d) and (e), or amalgams thereof. 

j. System B strictly noninter acting: Here system B 
is strictly noninteracting. Diagrams of the form (d) and 
(e) are not present in this scenario, and consistent with 
this the degenerate node survives. 

k. Hubbard interactions with weak intersystem hop- 
ping: Here the Green's function for system A + B is 
calculated with Hubbard interactions U nm = S nm U and 
V = t/3 rather than V = t as in the other cases. This 
scenario has the weakest overall coupling between the 
systems, and consistent with this, the smallest splitting 
of the degenerate node. However, the node splitting is 
still significant. 

Figure [5| A degenerate node from fine-tuned inter- 
actions: Finally, we present an example of a degener- 
ate node that exists within an interacting system despite 
the presence of diagrams that are inconsistent with se- 
ries propagation. The cost is that instead of its existence 
being ensured by the geometry of the system, the de- 
generate node requires fine-tuning of, for example, an 
interaction matrix element. Moreover, this example also 
demonstrates that a perturbation can destroy a degener- 
ate node by of lifting rather than splitting. 

We consider the system A + B with the parameter- 
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FIG. 10. The absolute matrix elements of th e family of single- 
particle Hamiltonians defined by Eq. (A3). The parameters 



of the four-site extended Hubbard model are tuned such that 
if 14 « near E — 0. Since H11, if 22, and H13 are sup- 
pressed at this point by particle-hole symmetry, the system 
has a "low-energy topology" that exhibits a degenerate node 
in the element G23 of the associated Green's function. This 
approximation breaks down near nodes in the eigenvalues of 
the Green's function (singularities in the self-energy), as can 
be seen around ±8 eV. 



ization given by U = U a = Up = 5 eV, U 1 = eV. 
In this case the supernode present within noninteract- 
ing models is not split, but completely lifted. This per- 
sists as C/14 = U 1 is increased, up to a critical value of 
w 1 eV near which there is a sharp transition between 
the regimes of node lifting and splitting. At the critical 
value a supernode exists despite the presence of interac- 
tions that are inconsistent with series propagation, but 
at the cost of fine-tuning. 

To investigate this phenomenon, it is instructive to use 
the exact Coulomb self-energy Y> c (E) to define a formal 
mapping of the many-body system A + B onto a family 
of single-particle systems with Hamiltonians: 



H 



E c (£) 



(A3) 



where H^\_ B is the noninteracting portion of the Hamil- 
tonian of system A + B. This reproduces the exact 
one-body Green's function by construction. The ele- 
ments of H are plotted in Figure (10) for the critical 



value U 1 = 0.969 eV. Regardless of tuning, elements cor- 
responding to next-next nearest neighbor hopping (e.g. 
Ef 3 ) break particle-hole symmetry and are therefore sup- 
pressed at half-filling near the particle-hole symmetric 
point (E = 0). At the critical value of J7 7 , nonlo- 
cal exchange cancels exactly with higher order diagrams 
and the element E^) also vanishes, whereas for larger 
(smaller) values of U 1 it is respectively of the same or 

opposite sign as the hopping matrix elements in H^]_ B - 
Thus, in this case, the competition between diagrams 
(d) and (e) leads to the switchover between node lifting 
and splitting depicted in Figure |5j Conclusive proof of 
the important qualitative role of higher order processes 
is furnished by self-consistent Hartree-Fock calculations 
(not shown), which do not lead to node lifting for this 
model. 
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FIG. 11. The element G23 of the Green's function for system 
A + B evaluated at the particle- hole symmetric point, plotted 
as a function of the average interaction matrix element U and 
a hopping matrix element A connecting sites one and four. 
For simplicity, we have taken V — t and constant interactions, 
i.e. Unm = U. A degenerate node exists for a family of 
parameterizations corresponding to the arc-shaped path that 
intersects U — eV, A = eV. However, it is forbidden 
by symmetry in the special case of A = t. Absent the exact 
symmetry present in this case, the interplay between topology 
and interactions is complex. 



As an aside, we note here that the approximate Hamil- 
tonian obtained by setting E = in the exact expression 
( |A3| ) is valid only away from nodes in the eigenvalues of 
the Green's functions (singularities in the Coulomb self- 
energy) . 

Figure\Tl\ A degenerate node from fine-tuning between 
a hopping matrix element and a charging energy U with 
constant interactions: The generic dependence of de- 
generate nodes upon fine-tuning in the presence of inter- 
actions can be seen by considering the amplitude G2z(E) 
evaluated at E = as a function of parameters in the 
Hamiltonian. For simplicity, this is done in Figure [TT] for 
system A + B with all the interaction matrix elements 
set to a constant U, with V = t, and a with a hopping 
matrix element A introduced between sites one and four. 
A degenerate supernode traces a path that corresponds 
to a family of models wherein H14 = 0. Just below this 
region the node is split, and above it, lifted. 

The case A = t is also special in that symmetry then 
precludes the existence of a node at E = regardless 
of the charging energy U. This can be argued formally 
based on the eigenvectors of the Green's function, which 
here are just equal to the noninteracting ones. More gen- 
erally, the same claim holds for non-constant interactions 
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FIG. 12. The absolute elements of the retarded Green's func- 
tion of a variation on the extended Hubbard system A + B, 
calculated at half- filling via exact diagonalization. Here the 
on-site energy of site 3 in system B is shifted up by 3 eV so 
that systems A and B are not identical and, moreover, there 
is no particle-hole symmetry. Panel (a) shows the results for 
pure series propagation (no diagrams from Figure [3]) , panel 
(b) shows the results with the inclusion of diagrams (a)-(c), 
and panel (c) shows the exact solution. 



because in the presence of such symmetry, the interacting 
eigenvectors can be taken as equal to the noninteracting 
ones. 

Figure[T^ Broken symmetry between A and B: Up to 
this point we have only considered cases wherein systems 
A and B are identical. If this is not so, then it should 
perhaps not be suprising that a degenerate node expected 
on the basis of the alignment of nodes in A and B is 



disrupted. Figure [12] explores this, depicting the Green's 
function when the on-site potential of site 3 is shifted up 
by 3 eV: 

Here panel (a) has only diagrams of the form (c), 
which, as before, do not perturb the node structure in 
qA,b ^ na k ghown) and therefore do not disrupt the nodes 
in G A+B . Panel (b) has all diagrams of the form (a)-(c) 
and collectively these do split the degenerate node now. 
Since series propagation is enforced here, this must occur 
via a disruption of the relevant nodes in G A,B . Moreover, 
to cause splitting these nodes must be shifted by different 
amounts, which requires that A and B be distinct. In this 
case it also necessitates broken particle-hole symmetry so 
that the nodes are not pinned at E = 0. 

It is interesting to note, however, that diagrams (a)-(c) 
still only make a small contribution to the overall node 
splitting, as can be seen in panel (c), which shows the 
exact solution. Thus, at least in this case, the disruption 
of the degenerate node is not due primarily to changes 
in the electronic structure of A or B as denned by the 
dressing of G A,B . 



3. Greater than four sites 

We discuss only briefly cases involving more than four 
sites. In particular, we consider an eight-site extended 



Hubbard model depicted in Figure [13] together with 
its node structure. Here the positions of local minima 
(dashed lines) and nodes (solid lines) in the element G27 
of the retarded Green's function are shown as a function 
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of a dimensionless prefactor 7 that premultiplies the in- 
teraction matrix U nm . When this prefactor is unity, the 
parameters in the Hamiltonian are consitent with those 
for an organic molecule with the appropriate symmetry. 

With 7 = this system exhibits a fourth order degen- 
erate node where G27 oc E A . This splits into two local 
minima as interactions are introduced. As the strength 
of interactions passes some critical value j c these minima 
bifurcate and produce two lowest order nodes. Each of 
these three cases (i.e. 7 = 0, < 7 < j C: and 7 > j c ) 
can be understood in terms of the properties of the roots 
of the polynomial from Eq. (10). 
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FIG. 13. The location of the nodes (solid lines) and minima 
(dashed lines) in the element C27 of the retarded Green's func- 
tion associated with the depicted eight-site extended Hubbard 
system. Here 7 is a dimensionless parameter that scales the 
strength of the interactions. The calculation was performed 
via exact diagonalization. 



t < 0. In particular, although the intuitive case 
of Gp a (t,0) = —ihQ(t)(d^(t)dl ( (0)) was considered 
first to formulate a definition of series propaga- 
tion, the preceding applies equally to the retarded 
Green's function 3 - 4 ^ from the Keldysh formalism: 
G a p(t,0) = -ihG(t)({dp(t),dl(0)}). Transport related 
quantities are elegantly expressed via this quantity, as it 
combines the amplitudes for particle-like and hole-like 
processes, both of which contribute to e.g. charge and 
heat transport. 

For a nanostructure wherein orbitals /i and v are cou- 
pled to a metallic electrode via hopping matrix elem ents, 
the corresponding tunneling self-energy is of the forrr i 35 * 43 -* 
Ej m = YfiSpnSum + T^8 vn 5 vm . Under these conditions 
the Dyson equation for the dressed retarded Green's func- 
tion: ePEr 



G = G + GTj T Q 



implies: 



Guv = 



(l-G^l)(l-G^ v ) 



(Bl) 

Here we have assumed that the full self-energy asso- 
ciated with coupling is equal to the tunneling self- 
energy, i.e. take n the so-called elastic cotunneling 
approximationP^MIH] 



In the broad-band limrtP^ Y/F 



and ^ 



in terms of elements and T v of a tunneling-width ma- 
trix. In this case the transmission function be tween the 
electrodes coupled to /i and v is then given by! 43|52 t ^ 



Appendix B: The connection between propagation 
and transmission 

In section (TT] we assumed little about the 
Green's function except that it vanishes when 



7V = FT"|a^| 2 ocG^ 

With regard to this work, the most important implication 
of the foregoing is that a node in G^ v is sufficient for a 
node in G^v and consequently T^ v . 



Appendix C: A two-pole approximation for the eigenvalues of the retarded Green's function in extended 

Hubbard models 



The retarded Green's function of an intera cting system can be expressed in terms many-body energies and eigen- 
states via a Lehmann representation. ^ 5 * 55 * 56 * In particular, at zero temperature: 

r = j_ v <*KI* / )<* , I'4I*) | 1 v (MaM rri x 

nm n N ^ E- [En,, - + i0+ n N E+ [E*, - E*] + i0+ { ' 

*'eB N+1 *'eB N _ 1 



Particle-like propagation Hole-like propagation 
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where B% = {|^)} is an orthonormal set of statistically occupied TV particle ground states, Bn±i = is an energy 

eigenbasis for the space of N ± 1 particle states, and and E^> are many-body energies. In the present context 
d) n creates a particle on the nth spin-orbital in an extended Hubbard model, and Qn is the number of statistically 
accessible N particle ground states. For the sake of brevity, we consider the spin degrees of freedom to be implicit in 
the index n. 

As an alternative to working in a localized basis, one can consider creation operators d) u {E) defined by: 

4 = $>H4(£) (C2) 

or equivalent ly: 

dt(E) = 5>l»>4 (C3) 

Here {\v(E))} is an energy-dependent single-particle basis chosen to diagonalize the retarded Green's function, i.e. 
such that: 

G = ^A»H (C4) 

V 

where \ V (E) are the eigenvalues of G. This is always possible for an isolated system because the retarded Green's 
function is a normal matrix^ in this case. Without loss of generality, we also assume that {\v(E))} are eigenvectors 
of a set of generators for all the one-body symmetries of the system under consideration. 



Formally, Eq. (CI) and the definition of {|^)} implies: 



1 g^—t^' \ cr^'— ^ 



Xv Sl N ^ n E- [Ey - E*] + *0+ + Q N ^ n E + [Ey, - E*] + i0+ ^° 5 ^ 



N—l 



Particle-like propagation Hole-like propagation 

where — \ I ^'\S V \^)\ 1 are spectral weights associated with the transitions between many-body states. Since 
all ^ G B% are associated with the same degenerate ground-state energy, this is equivalent to: 

x -= E E -si + io+ + T E £+7fW (C6) 



where: 

^ = £ £ «»> 

and: 

= .E"^/ — .E"^ (C8) 

In the sum above, when ^' corresponds to a correlation-induced excited state, it contributes a narrow resonance 
at high-energy, e.g. in the shaded region inaccessible at the Fermi level in Figures [6] and [8j In the white regions in 
the same figures, these can be neglected or, as seen shortly, treated by renormalizing spectral weights. This suggests 
that, in the absence of orbital degeneracy, the eigenvalues \ v be approximated as: 

Z! ' Z " (C9) 



" E-el + iO+ E + el + iO+ 
where: 

Z * = ^T E ( C1 °) 
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FIG. 14. The retarded Green's function of the four-site extended Hubbard model considered in section [A] calculated exactly or 
using zeroth or first order approximations for the eigenvalues of the Coulomb self-energy. The eigenvectors are evaluated with 
the energy fixed at E — 0. The linear approximation is essentially exact below E w ±8 eV, where the Coulomb self-energy 
becomes non-analytic due to a singularity associated with a node in an eigenvalue of G. 



and: 

z v = cT E ( cn ) 



Here, e^ ,h are energies associated with the particle-like and hole-like transitions with nonzero spectral weight closest 
to the Fermi level. Due to our choice of |z/), the foregoing is also valid in the presence of an orbital degeneracy if it is 
due to a symmetry not broken by interactions. 



To explore the approximation described by Eq. (C9), it is useful to consider briefly the limit of constant 
interactions, 58 i.e. U nm = U , in which case it is equivalent to its exact counterpart Eq. (C6). In this case the 



many-body states are Slater determinants and the eigenvectors \v) are just the noninteracting ones. The spectral 
weights are then zero or unity depending the occupancy of the spin-orbital \v) in |\J>) and \^ f ). Likewise, the 

energies e^ ,h are given exacly by = e v + U([N — No] + \) and e\ = — e v — U([N — No] — ~). There are therefore 



no correlation induced resonances at high-energies, and the resonances in Eq. (|C9j) for \ v correspond to the addition 
or removal of a particle from the single-particle state \v). 



If distance-dependent interactions are dialed on, i.e. U nm = U — > U + A/7 nm , the eigenvalues \ v in Eq. (C6) and 
potentially the eigenvectors {\v)} deviate from those with constant interactions. At low energies (e.g. the white region 
in Figures [6] and [8| , this comes primarily from a shift in the positions of the dominant resonances, which correspond 
to those present with constant interactions, as well as a reduction of their spectral weights Z^ h . The latter occurs as 
spectral weight is transferred to the correlation-induced resonances at high energy. 



Overall, we find that in the simple cases considered herein, an expression of the form ( C9 ) with "effective" parameters 
^P'* 1 , < ZP' h < 1, and U > reproduces the low-energy Green's function with high accuracy. Moreover, with long- 
ranged interactions these parameters are typi cally near to their bare values. While this method is similar to two-pole 
approximations applied to the Hubbard modeP^, it is expected to work best when the length scale of the interactions 
is larger than or comparable to the size of the system considered. 

This approxi matio n may also be cast in terms of the eigenvalues of the Coulomb self-energy E^, which, from the 
Dyson equatiorP^MI^ are related to the eigenvalues above simply by: 

A, = ^ (C12) 



where e v are the noninteracting energies. Eq. (C9) then corresponds to a Laurent expansion of the Coulomb self- 
energy with the form: 

£„ = c_i (E - Eo)' 1 + c + ci (E - Eo) (C13) 

If there is particle-hole symmetry so that = = \Z V , then Eo = e + U[N — Nq], c i = Z~ x •> 

Co = Eq — £ V1 and c\ = 1 — Z~ x . Here A^n is the number of sites in the extended Hubbard Hamiltonian. The residue 
is associated with the nodes in \ v and the linear term adjusts the spectral weight of the dominant resonances in 
in a manner reminiscent of a Fermi liquid PP. 

The assumption of particle-hole symmetry above is appropriate when \v) is expected to be half-filled. If, for example, 
\v) is instead expected to be completely empty, then Z^ = 0. In this case Eq = e v + U([N — No] + |), c_i = 0, 



17 



Co = Eq — s u , and c\ = 1 — Z" 1 . In general, this approximation breaks down near correlation-induced nodes in the 
eigenvalue, which correspond to singularities in the exact Coulomb self-energy that render it non- analytic. 

As a concrete example of the foregoing, we offer Figure [14} which shows the element G23 of the Green's function 
of system A + B (Figure [7]) with long-ranged Coulomb interactions as specified by Eq. (A2). This is calculated 
at half-filling using zeroth or first order approximations for the eigenvalues Tj v and with the eigenvectors \v) held 
constant. The linear (n = 1) approximation is excellent below E « ±8 eV, where the Coulomb self-energy becomes 
non- analytic due to a correlation-induced singularity in one of its eigenvalues. Equivalently, at this location there is 
a node in an eigenvalue of the Green's function. 
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